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We propose and study a "gold-washing" - type of algorithm which smooths out the short range fluctuations but 
leaves invariant instantons above a certain size. The algorithm needs no monitoring or calibration. 
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1. THE COOLING PROBLEM 

As a method to eliminate UV noise and permit 
studying topological excitations cooling must ful- 
fill certain requirements: 

- smooth out the short range fluctuations, includ- 
ing "dislocations" 

- preserve the structure at the physical scales, in- 
cluding size and location of instantons 

- ensure stability of the cooled structures 

- need no monitoring or any engineering which 
could make it configuration dependent or would 
only slow it down without ensuring stability. 

Since Wilson action has no stable instantons 
(they shrink and decay under cooling) one usually 
attempts to calibrate or engineer the cooling pro- 
cedure based on it such as to obtain metastability, 
after the noise has been reduced but before the 
instantons shrank to zero. This, however, makes 
cooling more an art than a method and also does 
not answer the question of preserving the physical 
scales of the original configuration. A general tree 
level analysis of various actions Q distinguishes 
between "under - improved" (e.g., Wilson) and 
"over - improved" ones; under the former instan- 
tons shrink, under the latter instantons beyond 
a certain size grow 0,^. Wc here consider the 
problem of "improved" cooling algorithms which 
ensure a high degree of scale invariance for the 
instantons beyond the short range scale. They 



should be useful not only for producing good sus- 
ceptibility data, but also for studying various con- 
jectured effects of instantons for the spectrum or 
the chiral transition by providing smooth config- 
urations which preserve the large scale structure 
of the original (hot) ones. This is why we called 
this approach "gold washing" . 

The question of a good cooling algorithm is of 
course related to that of an action possessing scale 
invariant classical solutions. "Fix point perfect 
actions", expected to have generally good scal- 
ing properties have been constructed in terms of 
many loops and higher representations One 
still needs, however, to test their usefulness for 
cooling algorithms in gauge theories. 

Our approach here has been different. Since 
we are primarily interested in instantons our 
construction refers directly to them. We start 
with tree level improvement and tune the ac- 
tion to obtain stable instanton solutions which 
are practically scale invariant beyond some small- 
size threshold; the latter should be such that only 
physically relevant instantons are preserved while 
short distance topological fluctuations ("disloca- 
tions" ) are eliminated in the cooling. It turns out 
that the improved cooling algorithm which will be 
presented below fulfills most of the requirements 
stated above. The instantons stabilized by it are 
well fitted by the continuum ansatz for volume in- 
tegrated quantities (we use twisted b.c. in "time" 
but the same should hold if the time extension 
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of the lattice is large enough). They have inte- 
ger topological charge within 0(0.1%). Higher 
charge configurations are well behaved under this 
algorithm, however instanton-antiinstanton (I-A) 
pairs annihilate of course after a certain num- 
ber of steps, due to their interaction. The study 
of pairs necessitates therefore further develop- 
ments. Although we concentrate here on the cool- 
ing problem, the good scaHng properties of the 
actions we study make them useful also as im- 
proved Monte Carlo actions. While our analysis 
here refers explicitely to SU{2), most relations 
are valid for SU (N) generally. 

2. IMPROVED COOLING 

One can take various loops in the action and 
tune their couplings to improve the approach to 
continuum. For simplicity we work only with fun- 
damental, planar loops of size m x n: 



minimizes the above local action at each step and 
involves no further calibration or engineering. 
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Figure 1. Action (a) and size (b) vs sweep number 
for various algorithms. S(W) is Wilson action. 
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= (l,l),(2,2),(l,2),(l,3),(3,3)for 
The Gibbs factor is exp {—-^S). 
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Here (m^, rii) 
i = 1,...,5. 
The choice 

5(e): ci = (4-e)/3, C2 = (e - l)/3 (3) 
C3,4,5 — 0, leads to "under-improved" actions for 
e > (instantons shrink, as for the Wilson ac- 
tion e — 1). For e < the actions are "over- 
improved" (instantons grow) and for e = there 
are no O(a^) corrections Using the results in 
Q one can construct a one parameter set of ac- 
tions that have no 0{a^) and 0{a'^) corrections 
ci = (19 - 55 C5)/9, C2 = (1 - 64 c5)/9 

C3 = (-64-^640 C5)/45, C4 = 1/5 - 2 cg (4) 
The simplest ones are S{ALi) (cs — 0) and 
S'(3I/i)(c5 = 1/10). Since for the configurations 
we have studied numerically the former tends to 
under-improve while the latter over-improves, we 
chose for our simulations a combination of both, 
S{5Li), with C5 = 1/20. This action showed al- 
ready good properties and we did not try to op- 
timize it further. The cooling algorithm exactly 



We first check the behaviour of various instan- 
tons under the improved cooling taking as start- 
ing configurations those obtained by cooling with 
5(e) with various e. All results are obtained on 
12"^ lattices and for gauge group SU{2). Since 
we wanted to disentangle small distance from fi- 
nite lattice size effects, we use throughout twisted 
p.b.c. with k — (1, 1, 1) twist, to ensure that in- 
stanton solutions exist on finite lattices |^J^ . Fig. 
1 compares the cooling behaviour of some algo- 
rithms {S(p) denote the underrelaxed algorithms 
of 1^): (a) the action (notice the different scales) 
and (b) the size Ppeak = 6'n'^Ns/S{t)\max for the 
same instanton, as function of the cooling step 
{Ns = 12 here). Fig. 2 shows the behaviour of 
instantons of various initial sizes under the 5Li- 
cooling. As it can be seen, they remain practically 
unchanged over any practicable number of cool- 
ing sweeps and for p > po — 2.3 (in units of a), 
which seems to represent a stability threshold for 
this cooling algorithm. The size of the original 
instanton has been varied by applying 5(e) cool- 
ing, therefore these instantons do not correspond 
to minima of 5(5Li). Submitting such configura- 
tions to 5Li cooling first readapts them to the new 
equations of motion, implying the small changes 
which can be observed in the first ~ 50 sweeps 
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Figure 2. Instanton size under 5Li improved cool- 
ing starting from various configurations. 
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Figure 3. Instanton evolution (action and size) 
during 5Li cooling. 



of Fig. 2. Fig. 3 shows the effect of 5Li cooling 
on instantons of different sizes taken from vari- 
ous stages of Wilson cooling of 3 different starting 
configurations (the upper curves give the values of 
S{5Li) for these Wilson cooled instantons). We 
see how the S{bLi)min emerges as an envelope. 
We observe again the stability threshold around 
2.3 and the extent of configuration dependence at 
small p. The triangles are data from Monte Carlo 
instantons (/3 = 2.5, 12'' lattice) cooled with 300 
sweeps of 5Li cooling. 

To describe the instantons we use the volume in- 
tegrated action and charge densities ( "profiles" ) : 
S{t) = NsJ2six,y,z,t) (5) 

x,y.z 

Q{t) = 8TT^NsyQ{x,y,z,t) 



x,y,z 



(6) 



to be fitted by the continuum ansatz with peri- 
odicity satellites, e.g. for the charge: 
Q{t) = ±6TT^N,p-'[f{t) + fit + Nt) + fit - Nt)] 

f{t)^[it-tof+pT'^' (7) 
Here to defines the location and p the width, 
"Pprofiie"- The latter agrees rather well with 
ppeak as long as the total charge is ±1. This is en- 
sured at not too small p, i.e. for stable instantons. 

Since Wilson action produces short range fluctu- 
ations one can ask whether the effective cut-off at 
Pq ~ 2.3 ensures the absence of unphysical fluc- 



tuations in cooled MC configurations. Following 
an argument of Pugh and Teper we write the 
contribution of small instantons as (iV = 2 here) 
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We can calculate Swipo) by cooling large instan- 
tons with the Wilson action (they shrink). One 
can then flnd yOo(a(/3)) from the condition that the 
two factors in Eq.(8) compensate each other such 
that the contribution from instantons larger than 
Po stays finite in continuum. We obtain a very 
flat dependence: po = 1.9276, 1.9280, 1.9284 for 
/3 = 3, 6, oo, respectively. Hence a threshold po of 
about 2 and independent on f3 seems satisfactory. 

3. MONTE CARLO ANALYSIS 

To illustrate the capability of improved cooling 
we apply 5Li cooling to a Monte Carlo simula- 
tion at /3 = 2.5 on a 12^ lattice, k = (1, 1, 1) twist. 
We use 160 configurations generated by heat bath 
with Wilson action and taken 250 sweeps apart 
(after 20000 thermalization sweeps) . In Fig. 4 we 
show results from 300 cooling sweeps. The sta- 
ble situation typically sets in between 50 and 100 
sweeps. Instantons very near below the thresh- 
old Po — 2.3 may need more sweeps to disappear 
- see Fig. 2; this introduces a systematic error 
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Figure 4. MC analysis: typical charge 1 and 3 
configurations and size and charge distribution. 
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Figure 5. Early evolution of a MC conf.. The line 
for 1 sw is an eye-guide, the rest fits based on (7). 



in the number of stable instantons in this region, 
which vanishes with increasing number of sweeps. 
In Fig. 4a,b we show two typical configurations, 
the first with one wide {pprofUe = 4.66) instan- 
ton and the second with total charge 3 showing 3 
instantons of widths 2.80, 4.47 and 3.37 located 
at t=1.22, 6.22 and 9.29, together with the fit 
Eq. (7) . In Fig. 4c wc show the size and in 4d the 
charge distributions. The latter is representative 
for physical instantons, i.e. beyond the thresh- 
old po- A gaussian fit cxp(— 6(5^) gives b = 0.34 
at cooling sweep 20 and b = 0.37 beyond 100. 
The charge stabilizes to within 1% an integer be- 
tween 20 and 100 sweeps, therefore susceptibility 
and charge distribution can be estimated quite 
early. The topological susceptibility extracted 
from data at 300 sweeps is = 6.8(9) 10"^. 

The size distribution, however, represents only 
the stable situation and does not take into ac- 
count TA pairs which would have annihilated be- 
fore. For illustration we show in Fig. 5 the charge 
and action history (a) and the evolution of a typi- 
cal configuration imder the 5Li cooling (b,c). The 
initial structure seems to have a double instanton 
and 2 I- A pairs. An I- A pair can still be seen with 
the continuum ansatz Eq. (7) at sweep 5 (sizes 
1.38 and 2.95) and at sweep 10 (sizes 2.74 and 
3.36), together with the double instanton (size 2.5 
at sweep 5 and 2.86 at 10) - see Fig. 5b. After 



sweep 20 (Fig. 5c) only the latter remains (size: 
3.07 at 20, 3.29 at 50, 3.52 at 300 sweeps). Hence 
this annihilation may result in a depletion of the 
lower part of the size distribution. The question 
of describing complex configurations is a separate 
problem currently under study. 
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